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En este artículo se presenta una revisión del conocimiento actual sobre las propiedades elásticas 
de las membranas de origen biológico. Bajo la hipótesis del mosaico fluido se considera a una bi- 
capa de fosfolípido como la base estructural de una membrana biológica y se expone un modelo 
debido a Helfrich a partir del cuál se obtiene una expresión para la energía elástica de una mem¬ 
brana. La energía libre queda completamente caracterizada por sus curvaturas locales principales y 
cuatro parámetros dependientes de la composición química de la membrana y de su entorno local. 
Asimismo, se expone la justificación teórica de un método experimental (que hace uso de microma- 
nipulación con micropipetas) que permite medir la constante elástica de curvatura media de una 
membrana. Finalmente se describen algunas consecuencias físicas de estos conceptos, como la ex¬ 
plicación de la transición de una fase lamelar a una esponja, o la aparición de la interacción estérica 
entre membranas. 

Descriptores: Energía libre de Helfrich; membranas biológicas; bicapas lipídicas; elasticidad de 
membranas. 

In this article we review the present knowledge on the elastic properties of membranes of biological 
origin. Assuming the fluid mosaic hypothesis, we consider the phospholipid bilayer as the structural 
base of a biological membrane. We expose a model due to Helfrich which can be used to obtain an 
expression for the membrane elastic energy. The free energy is completely characterized by two local 
principal curvatures and four parameters that depend on the Chemical composition of the membrane 
and its local environment. In addition we present the theoretical justification of an experimental 
method (the Micropipette Manipulation technique) that may be used to measure the membrane 
rigidity (the elastic constant that corresponds to the mean curvatura of the membrane). Finally, 
we describe some physical consequences of these concepts, such as the explanation of the lamellar- 
sponge phase transition and the emergence of the steric interaction between membranes. 
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I. INTRODUCCIÓN 

Las membranas biológicas son componentes primor¬ 
diales de todos los organismos di- Se trata de es¬ 
tructuras esencialmente formadas por una bicapa de fos¬ 
folípido con proteínas incorporadas, que están presentes 
en las células biológicas así como en los diferentes organe- 
los. 

Las membranas tienen diferentes funciones biológicas. 
Por una parte, aíslan a la célula misma y a ciertos or- 
ganelos, como el núcleo celular, de su exterior. Pero 
además de servir como fronteras entre los espacios in- 
tra y extracelular o inter-extra organelo, las membranas 


también participan en otras funciones como el transporte 
activo, el encapsulamiento de sustancias, el soporte de 
canales iónicos y la estabilidad de la propia célula. Las 
membranas también juegan un papel esencial en proce¬ 
sos biológicos como la endocitosis, la exocitosis, la fusión 
celular y la formación de protuberancias [í|, Q . 

Análisis químicos y estudios biológicos han permitido 
conocer la composición química, así como ciertas fun¬ 
ciones biológicas de las membranas. Sin embargo, ex¬ 
isten aún preguntas sin resolver en este campo de la 
Ciencia. Por ejemplo: ¿Cuáles son las propiedades 
físicas relevantes de las membranas biológicas? ¿Cómo 
dichas propiedades determinan o influencian la estruc- 
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Figura 1. Forma que los glóbulos rojos adoptan en su estado 
natural. Note como la membrana que rodea a cada célula se 
tiene que curvar para adecuarse a la forma celular. 


tura y la función biológica de estos agregados? En¬ 
tre las propiedades físicas que reciben especial atención 
en estos tiempos se encuentran las dinámicas, estruc¬ 
turales y mecánicas. En particular, entre las propiedades 
mecánicas destaca la elasticidad de las membranas, 
fenómeno que está relacionado de una manera apenas 
parcialmente conocida por el momento con fenómenos 
como la endocitosis y exocitosis, y que juega un papel 
relevante en la forma que adquieren algunos organelos en 
la célula (retículo endoplásmico, complejo de Golgi). 

En este trabajo, nos proponemos comentar algunos 
conceptos importantes para entender ciertas propiedades 
elásticas de las membranas. Prestaremos especial aten¬ 
ción a la teoría de Wolfgang Helfrich, científico alemán 
que llegó a la expresión utilizada actualmente para la 
energía elástica de una membrana. 

Históricamente, el estudio de las propiedades elásticas 
de membranas está relacionado con el problema de pre¬ 
decir teóricamente la forma de los glóbulos rojos (Figura 
[T|). Estas células, y por tanto la membrana que las rodea, 
tienen forma discoidal, con una concavidad característica 
en sus caras. En la búsqueda de entender el origen de esta 
estructura, en 1970 Canham @ propuso un modelo en el 
que se considera constante tanto el área como el espesor 
de la membrana celular. En dicho modelo, los glóbulos 
rojos están caracterizados únicamente por su elasticidad 
respecto a la curvatura, es decir, mediante una energía 
elástica de curvatura que, por unidad de área, adquiere 


la siguiente forma: 

/ = ( 1 ) 

donde k representaba el módulo elástico de curvatura me¬ 
dia de la membrana y H la curvatura media de ésta (am¬ 
bas cantidades serán definidas en la sección III Al) . Ver¬ 
emos más adelante que esta expresión es parcialmente 
correcta. Sin embargo, este modelo conduce a formas no 
observadas en los glóbulos rojos 0. 

En la misma década de los setentas, la comunidad de 
científicos de la materia condensada aprovechó las teorías 
relativamente maduras de la física de cristales líquidos 
termotrópicos para aplicarlas al estudio de membranas. 
Esto es posible debido a que localmente una membrana 
se comporta como una celda de cristal líquido nemático 
en alineamiento homeotrópico; esto es, moléculas con 
forma de “rodillos” alineadas perpendicularmente al sus¬ 
trato. En este tipo de materiales, las moléculas tienen 
un orden orientacional de largo alcance. 

Aprovechando estos hechos, Helfrich utilizó la ex¬ 
presión para la densidad de energía libre de Frank para 
cristales líquidos uniaxiales considerando la normal 
a la membrana como el director del cristal líquido. Así, 
dedujo la energía de curvatura por unidad de área de 
una membrana @ . La ecuación resultante, publicada en 
1973, es llamada "energía libre de Helfrich” para mem¬ 
branas fluidas y es la base del análisis de las propiedades 
elásticas de las membranas. En este trabajo presentare¬ 
mos una deducción alternativa de la energía de Helfrich 
y comentaremos algunas de sus consecuencias físicas. 

Nuestro artículo está dividido de la siguiente manera. 
En la sección \LK\ se comentará brevemente la ubicación 
biológica, así como la estructura química de las mem¬ 
branas celulares. En la sección El se expone el modelo 
de Helfrich para biomembranas. En las secciones sigu¬ 
ientes se exponen algunas de las consecuencias físicas 
de la elasticidad de curvatura: fluctuaciones térmicas 
de membranas (sección IH Ep . métodos de medición de 
las constantes elásticas (IHFI) . transformación de la fase 
lamelar a la fase esponja (IHI Al) e interacciones estéricas 
entre membranas (IHIBI) . Finalmente, en la sección final 
damos algunas conclusiones del trabajo. 


A. Agregados y bicapas lipídicas 

Las membranas biológicas son estructuras muy com¬ 
plejas formadas por una gran variedad de moléculas: 
fosfolípidos, proteínas, etc. [í|. Estudios basados 
en microscopía electrónica y en el análisis de la com¬ 
posición química de biomembranas, así como estudios 
de propiedades físicas tales como la permeabilidad y la 
difusión de proteínas y de moléculas lipídicas en mem¬ 
branas, condujeron al desarrollo del llamado modelo 
del mosaico fluido. Este modelo, propuesto por S. J. 
Singer y G. Nicolson en 1972 Q, es el más aceptado 
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para describir a las membranas biológicas. La estruc¬ 
tura base de la membrana es una bicapa lipídica donde 
las moléculas individuales pueden moverse como en un 
líquido bidimensional. La bicapa lipídica es una mezcla 
de varias clases de moléculas, en particular fosfolípidos 
y glicolípidos, más otras moléculas pequeñas como coles- 
terol. Las proteínas de membrana se encuentran incrus¬ 
tadas en la bicapa lipídica. La estructura de las mem¬ 
branas biológicas es posible debido a que los fosfolípidos 
pertenecen a un tipo especial de moléculas, llamadas “an- 
fifílicas”, a la cual también pertenecen los surfactantes. 
Este nombre se debe a que una parte de la molécula (la 
cabeza polar) es soluble en agua, mientras que la otra 
parte (la cola hidrofóbica) es insoluble en dicho solvente. 
Así, cuando estas moléculas se disuelven en agua, las 
partes hidrofóbicas tienden a asociarse para evitar el con¬ 
tacto con el solvente, donde son insolubles, lo cual da lu¬ 
gar a la estructura en forma de bicapas (recordemos que 
el medio celular es esencialmente acuoso). Las proteínas 
que se pueden insertar en dicha bicapa, tienen segmentos 
hidrofóbicos que tampoco pueden solubilizarse en agua, 
los cuales se asocian con el centro hidrofóbico de la mem¬ 
brana. En la Figurarse muestra de forma esquemática 
la estructura de una bicapa lipídica, así como la fórmula 
molecular de uno de los fosfolípidos más conocidos, de la 
familia de las fosfatidilcolinas. 

Mediante diferentes técnicas experimentales, como la 
dispersión de radiación y la microscopía electrónica 

[l| , se ha podido determinar el espesor de una membrana 
formada por moléculas anfifílicas (alrededor de 50Á). Di¬ 
cho espesor experimental coincide bastante bien con lo 
esperado a partir de la longitud de las moléculas que for¬ 
man las membranas. 

Es importante subrayar que algunas propiedades físicas 
(mecánicas y termodinámicas) de las membranas de fos- 
folípido dependen crucialmente del tipo de lípido. Esto 
se refleja principalmente en la denominada temperatura 
de transición, T^, de los fosfolípidos. Por debajo de esta 
temperatura, las moléculas del fosfolípido adquieren una 
estructura tipo gel (sólido amorfo bidimensional). Para 
temperaturas mayores que T^, la membrana es fluida. 
Experimentos de calorimetría han mostrado cómo de¬ 
pende de las características químicas de los lípidos m- 
Por ejemplo, para una cabeza polar dada (como fosfa- 
tidilcolina), Tm se incrementa conforme se incrementa el 
número de átomos de carbono en la cadena hidrofóbica. 
Esto se debe a que las cadenas de longitud grande magnif¬ 
ican el efecto hidrofóbico, dando como resultado mayor 
cohesión entre las moléculas que forman la membrana. 
De manera similar, la presencia de insaturaciones (dobles 
enlaces químicos) en las colas hidrofóbicas hace que Tm 
disminuya, efecto que depende también de la posición 
en que se presenta la insaturación m- Se ha observado 
también que la presencia de cierta asimetría entre las dos 
cadenas que forman la cola hidrofóbica de la molécula, 
hacen que Tm disminuya [lTj |. Esto tiene como conse¬ 
cuencia que los fosfolípidos más utilizados para preparar 
modelos experimentales de membranas biológicas sean 


aquellos con una cadena saturada y otra insaturada, 
como es el caso del l-estearoil-2-oleil-glicero 3-fosfocolina 
(SOPC), o de mezclas como la lecitina de soya o huevo. 
Esto asegura que la temperatura de transición, T^, esté 
siempre por debajo de la temperatura a la cual se llevan 
a cabo los experimentos. 

Cabe señalar que las moléculas anfifílicas no sola¬ 
mente son capaces de formar las membranas biológicas, 
sino también otro tipo de estructuras. Cuando estas 
moléculas (ya sean fosfolípidos o surfactantes) se disuel¬ 
ven en agua, su característica dual en cuanto a la solu¬ 
bilidad hace que se puedan formar diferentes estructuras, 
como las micelas (agregados individuales que pueden ser 
esféricos o cilindricos) o las bicapas. A su vez, las bicapas, 
que son la estructura básica de una membrana, pueden 
adoptar configuraciones en forma de vesículas (agregado 
esférico), fase lamelar (membranas paralelas), fase cúbica 
(membranas en arreglo periódico) o fases esponja (bica¬ 
pas desordenadas pero conectadas). Una representación 
esquemática de algunas de estas fases se muestra en la 
Figuráis] 

Cada molécula particular presenta algunas de estas 
fases en diferentes condiciones fisicoquímicas: concen¬ 
tración, temperatura, aditivos, etc. Por otra parte, se 
han desarrollado algunos modelos teóricos que predicen 
la topología de los agregados. El más sencillo de di¬ 
chos modelos, propuesto por Israelachvili [I^, estipula 
que la forma del agregado está determinada por el lla¬ 
mado “parámetro de empaquetamiento”, el cual se rela¬ 
ciona con la geometría de la molécula. Según este mod¬ 
elo simple, la asociación de moléculas de forma cónica da 
lugar a un agregado con gran curvatura, como las mice¬ 
las. En cambio, la asociación de moléculas relativamente 
cilindricas da lugar a agregados con poca curvatura, como 
las membranas. Este modelo permite explicar algunas 
transiciones observadas en las fases de moléculas an- 
flfílicas. Por ejemplo, si se tiene un surfactante cargado 
que forme micelas esféricas, esto es debido a la forma 
efectivamente cónica de la molécula, magnificada por la 
repulsión electrostática entre cabezas polares. Al agregar 
sal al sistema, las interacciones electrostáticas se apan- 
tallan, lo cual tiene el efecto de disminuir el tamaño efec¬ 
tivo lateral de las cabezas polares, por lo que la molécula 
se asemeja más a un cilindro que a un cono. En estas 
condiciones, la teoría del parámetro de empaquetamiento 
predice que la estructura que forman dichas moléculas 
evoluciona de la forma micelar hacia la forma de mem¬ 
brana con el sólo hecho de agregar sal al sistema. Este 
efecto ha sido observado experimentalmente. 

Sin embargo, hay otras transiciones que no han po¬ 
dido ser explicadas completamente mediante modelos 
geométricos tan sencillos. Una de ellas es la transición de 
la fase lamelar a la fase esponja (Figura [3]), donde mem¬ 
branas que tienen una configuración relativamente plana, 
adquieren una topología local en forma de silla de mon¬ 
tar a caballo. Se trata de una transformación topológica 
de las membranas. Lo mismo puede decirse en cambios 
donde se pasa de una fase lamelar a una fase vesicular, 
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Figura 2. Una de las componentes principales de una membrana biológica es una matriz lípidica, en la cual, las colas hidrofóbicas 
de las moléculas de fosfolípido evitan el contacto con el agua. En la figura mostramos también la molécula de uno de los 
fosfolípidos más comunes, de la familia de la fosfatidilcolina. 



Figura 3. Representación de la estructura de algunas de las 
fases típicas formadas por moléculas anfifflicas. La estructura 
más sencilla es la micela. Por su parte, las membranas pueden 
ordenarse de maneras diferentes, tal como en la fase lamelar, 
en la esponja o en la vesicular. 

mediante la adición de un polímero por ejemplo, lo cual 
también implica un cambio topológico. Veremos que este 
tipo de transformaciones se pueden entender en términos 
de la energía elástica de membranas, concepto que se dis¬ 
cutirá en la siguiente sección. 


II. MODELO DE HELFRICH 

Algunas propiedades geométricas y termodinámicas de 
las bicapas anfifílicas pueden ser entendidas desde un 
punto de vista unificado considerando la energía libre 


asociada a sus deformaciones. El espesor tan pequeño 
de una membrana (50Á), al menos comparado con las 
dimensiones de una célula (del orden de mieras), permite 
aproximar a la bicapa como una superficie bidimensional 
inmersa en el solvente M- 

La membrana se encuentra sujeta a las fluctuaciones 
debidas a la energía térmica del solvente. Si suponemos 
que la bicapa se encuentra en equilibrio con la solución 
acuosa, entonces, su energía está caracterizada por las 
deformaciones que la membrana puede sufrir; a) la com¬ 
presión o expansión de la bicapa en la dirección lateral y 
b) las deformaciones en la dirección normal a la superfi¬ 
cie (esto es, por la curvatura local de la superficie) [l3| . 
Bajo estas condiciones, la energía libre (por molécula) 
asociada a la membrana, puede ser expresada como una 
función del área promedio de cada molécula, E, y de las 
curvaturas locales asociadas a la superficie. 

Los conceptos matemáticos necesarios para llegar a la 
energía elástica de las membranas pertenecen a la Ge¬ 
ometría Diferencial de Superficies, por lo que conviene 
repasar algunos conceptos básicos de dicho campo de la 
matemática. 


A. Conceptos matemáticos previos 

Considérese un plano P, perpendicular al plano tan¬ 
gente (en un punto cualquiera) a la superficie que repre¬ 
senta a la membrana. A la intersección de la superficie 
con el plano P se le llama una “curva normal”. Esta 
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curva normal puede ser una línea recta, sin curvatura, 
o puede poseer cierta curvatura, a la cual se le llama 
“curvatura normal”, c. Estos conceptos se ilustran en 
la Figura SI Para asociar un valor a la curvatura nor¬ 
mal, notemos que una línea recta (por definición sin cur¬ 
vatura), puede verse como la circunferencia de un círculo 
de radio infinito. Esto sugiere definir a la curvatura como 
el inverso del radio de curvatura: c = 1/i?. Como es de 
esperarse, para una línea recta c —>■ 0 puesto que i? —>■ oo. 
A la curvatura así definida se le asocia un signo para dis¬ 
tinguir las dos posibles orientaciones de la curva: concavi¬ 
dad hacia un lado o concavidad hacia el otro, direcciones 
a las cuales se asigna arbitrariamente el signo positivo o 
negativo. 

Notemos que, dado que el plano P se definió simple¬ 
mente como perpendicular al plano tangente, existe una 
infinidad de posibles orientaciones para P. Dependiendo 
de dicha orientación, las diferentes curvas normales posi¬ 
bles tienen diferentes curvaturas. Sin embargo, puede 
demostrarse que hay dos direcciones particulares, per¬ 
pendiculares entre sí, para las cuales la curvatura normal 
es un extremo (en un caso máximo y en el otro mínimo). 
Dichas curvaturas normales extremas son llamadas las 
“curvaturas principales” de la superficie y se denotan me¬ 
diante los símbolos ci y C 2 . Todas las demás curvaturas 
normales en el mismo punto (así como la curvatura local 
de cualquier curva contenida en la superficie y que pase 
por dicho punto) pueden ser expresadas en términos de 
las dos curvaturas principales. En otras palabras, todas 
las deformaciones locales de la superficie en la dirección 
normal al plano tangente pueden ser caracterizadas a 
primer orden por las curvaturas principales ci y C 2 . Es¬ 
tos son los dos parámetros necesarios para caracterizar 
la curvatura de una membrana fluctuante. 

Sin embargo, es más conveniente matemáticamente re¬ 
alizar un cambio de variable, de ci y C 2 , a dos combi¬ 
naciones de ellas, la llamada “curvatura media”, P[ = 
(ci -I- 02 ) 12 , y la llamada “curvatura gaussiana”, K = 
C 1 C 2 . H es el parámetro mencionado al comentar el mod¬ 
elo de Canham en la introducción (ecuación 1). De he¬ 
cho, a segundo orden en las curvaturas principales, \H\ 
y K son las únicas cantidades invariantes ante cualquier 
parametrización de la superficie. Es importante señalar 
que, dado que tanto ci y 02 tienen asociado un signo, el 
valor de H depende de la orientación elegida como posi¬ 
tiva para la curvatura (cambiar la convención de signos, 
cambia el signo de ií). El signo de K, por su parte, es 
independiente de esta elección. La Figura [S] muestra al¬ 
gunas superficies con valores diferentes de las curvaturas 
media y gaussiana. 


B. Derivación de la energía de Helfrich 

Apliquemos ahora los conceptos geométricos antes 
mencionados a la descripción física de las bicapas lipí- 
dicas. La energía libre de la bicapa (por unidad de área), 
/, debe ser invariante bajo cualquier parametrización de 


la superficie, dado que es un concepto físico que no de¬ 
pende de la descripción matemática utilizada. Luego en¬ 
tonces, dicha energía debe ser una función de ií y K. 

En el equilibrio las bicapas pueden considerarse aprox¬ 
imadamente planas, al menos localmente; las deforma¬ 
ciones respecto a dicho estado están caracterizadas por 
curvaturas muy pequeñas; es decir, el radio de curvatura 
correspondiente es mucho mayor que el espesor de la bi¬ 
capa. Esto permite desarrollar la función / = f{H, K) 
mediante una serie de Taylor (a segundo orden en las 
curvaturas Ci y C 2 ): 

f{H, K)^fo + fiH + /aíf" + hK (2) 


donde las constantes /o, /i, ¡2 y /a dependen en prin¬ 
cipio del área promedio por molécula; en particular /o 
corresponde a la energía libre de una bicapa plana. 

Nótese que el término lineal en P[ depende de la ori¬ 
entación de la superficie. Cuando las monocapas que 
forman la bicapa tienen la misma composición, no hay 
razón alguna para dicha dependencia en la orientación y, 
en tal caso, la constante /i debe ser cero. En general, 
sin embargo, existe la posibilidad de una asimetría entre 
las monocapas y es precisamente el término lineal en P[ 
el que toma en cuenta este efecto. 

Completando el binomio cuadrado respecto a H, la 
ecuación © puede reescribirse como: 


Si entonces definimos las cantidades: 


CT = /o - 


« = 2/2, 


= _A 

2/2 

K = /s 


IL 

4/2’ 


la energía libre, por unidad de área, de la bicapa toma la 
forma 


f = ]^k{H - cq)^ + kK + a, (3) 

donde a, k, R y cq son la tensión superficial, el módulo 
de curvatura media, el módulo de curvatura gaussiana 
y la curvatura espontánea de la membrana, respectiva¬ 
mente. Colectivamente esas variables son conocidas como 
los “parámetros de Helfrich”. Dichos parámetros son 
independientes de la forma geométrica de la membrana 
pero dependen en principio de su composición química. 

La energía total de la bicapa (conocida como la energía 
de Helfrich) puede ser calculada integrando sobre el área 
total de la misma: 


F = 


1 9 

-k{P[ — co) + kK + a 


dA. 


( 4 ) 
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Figura 5. Superficies con distinta curvatura media y gaus- 
siana. 


Wolfgang Helfrich llegó a estas ecuaciones en 1973 en 
un afán de describir la elasticidad de las membranas 
biológicas La derivación original publicada ese año, 
se basa en una analogía con la energía libre para cristales 
líquidos propuesta por Frank Q. Sin embargo, la forma 
en que aquí arribamos a tal relación fue desarrollada por 
Safran en 1994. Hemos elegido esta última derivación 
debido a la sencillez y elegancia. 

Es importante notar que si no hay fuerzas externas que 
tensen la membrana, el área de la bicapa puede consider¬ 
arse constante. Esto es así porque las colas hidrofóbicas 
de las moléculas de fosfolípido son casi totalmente in¬ 
solubles en agua y están agrupadas de manera que el 
área por molécula se mantiene prácticamente constante. 
Modificar dicha área requiere de energías mayores a la 
necesaria para simplemente reacomodar las moléculas de 
la bicapa de modo que ésta se curve. Por esta razón, en 
ausencia de fuerzas externas, el término de la energía de 
Helfrich que contiene la tensión superficial, ct, es una con¬ 
stante sin relevancia al momento de minimizar la energía 
respecto a las posibles deformaciones de la superficie. En 
la discusión subsiguiente ignoraremos este término, ex¬ 
cepto cuando se mencione explícitamente lo contrario. 


C. Significado físico de los parámetros en la 
energía de Helfrich 

El significado de los parámetros de Helfrich puede en¬ 
tenderse sencillamente analizándolos por separado. 

Para comentar la física que se desprende del módulo 
de curvatura media, k, es conveniente empezar haciendo 
algunas observaciones acerca del término de la curvatura 
gaussiana, el segundo de la ecuación El llamado teo¬ 
rema de Gauss-Bonnet muestra que este término es 
una constante si la topología de la membrana no cambia. 
En efecto, este teorema estipula que para una superfi¬ 
cie cerrada de topología fija, la integral de la curvatura 
gaussiana es una constante, 

J KdA = Att (1 — p) , (5) 

donde p denota el género de la superficie, concepto 
topológico que está relacionado con la cantidad de poros o 
perforaciones que tiene la superficie. Algunos ejemplos de 
superficies con distinto género se muestran en la Figura 
[ü] La conclusión principal del teorema de Gauss-Bonnet 
es que, al estudiar la estructura de mínima energía de 
membranas de topología fija, el término que involucra 
la curvatura gaussiana puede ser ignorado, pues es con¬ 
stante mientras la topología no cambie. En otras pal¬ 
abras, la energía de curvatura asociada al módulo gaus- 
siano no cambia mientras la topología de la membrana 
sea la misma. 

Bajo las condiciones anteriores, la energía elástica de 
Helfrich toma la forma de un potencial armónico F = 

-k/ {H — co)^dA que depende de la curvatura media Fl 

en forma análoga a como la energía potencial de un re¬ 
sorte depende de su deformación. Entonces, el módulo de 
curvatura media, n, es análogo a la constante de rigidez 
de un resorte, y puede interpretarse como una constante 
de rigidez (ante deformaciones de curvatura) de la mem¬ 
brana. Guando k tiene valores pequeños, la membrana es 
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flexible; cuando k asume valores grandes, la membrana 
es rígida. Entonces, el término de curvatura media está 
relacionado con las fluctuaciones de forma u ondulaciones 
(sin cambio topológico) que experimenta la membrana 
por efectos de la energía térmica del solvente. 

Por su parte, la curvatura espontánea cq es la cur¬ 
vatura media que adquiere la bicapa cuando está libre 
de tensión, jugando un papel análogo al de la longitud 
de equilibrio en el caso del resorte. Hicimos notar an¬ 
teriormente que para monocapas idénticas la constante 
/i de la ecuación (¡2) debe ser cero y, por tanto, cq debe 
ser nula también. Es claro entonces que la curvatura 
espontánea está asociada a la asimetría de las monoca¬ 
pas que forman la bicapa. Esta curvatura es debida al 
hecho de que en general las monocapas se forman por 
distintas moléculas anfifílicas ocasionando que las bica- 
pas tiendan a curvarse espontáneamente hacia alguno de 
sus lados cuando dos de estas monocapas se superponen. 
Por supuesto que cuando la bicapa está formada por un 
solo tipo de moléculas, no existe curvatura espontánea y 
el estado de equilibrio tiende a ser el de una membrana 
plana. En resumen podemos decir que la ecuación ([3]) 
establece que cuando no hay cambios de topología en el 
sistema la curvatura media que minimiza la energía li¬ 
bre corresponde a la curvatura espontánea. El radio de 
equilibrio de la membrana es Rq = I/cq. 

Finalmente, para tener una idea del significado físico 
del módulo de curvatura gaussiana, k, imaginemos una 
bicapa formada por monocapas idénticas (cq = 0) que 
puede cambiar de topología. El término de la curvatura 

media, F = —k J H^dA, es mínimo para una superficie 

con H = 0, que puede ser un plano pero también una 
superficie localmente de tipo “silla de montar” (Ver figura 

0 . 

Puede ahora notarse que el término de la curvatura 
gaussiana, f RKdA, favorece superficies con distintas 
topologías, dependiendo del valor del módulo elástico 
gaussiano. En efecto, cuando k > 0, la energía rela¬ 
cionada con el término gaussiano se minimiza si los ra¬ 
dios de las curvaturas principales tienen diferente signo, 
es decir, si la curvatura gaussiana es negativa, como en 
el caso de superficies tipo “silla de montar”, prevale¬ 
cientes en las fases esponja y cúbicas. Por otra parte, 
cuando k < 0, el término gaussiano favorece superficies 
de curvatura gaussiana positiva (radios de curvatura en 
la misma dirección), contribuyendo a la estabilidad de 
agregados cerrados tales como vesículas. 

Esta breve discusión nos permite vislumbrar el tipo 
de fenómenos físicos que afectan a una membrana y que 
pueden ser descritos con la ecuación ©: forma de equi¬ 
librio, fluctuaciones térmicas y cambios de topología. Por 
esta razón, la energía de Helfrich ha sido utilizada en 
una gran cantidad de trabajos que describen algunas 
propiedads físicas de las membranas. De hecho, la posi¬ 
bilidad de calcular la energía libre de una bicapa ha per¬ 
mitido modelar sistemas biológicos reales como el aparato 
de Golgi [l^ o la forma discoidal de los glóbulos rojos 



Figura 6. Superficies con distinto género. 

i El- La forma en que se usa la relación m en los 
distintos modelos consiste en sumar las diferentes con¬ 
tribuciones energéticas para después minimizar la energía 
total a fin de obtener información sobre la forma o las 
fluctuaciones del sistema. En la siguiente sección vere¬ 
mos una forma matemáticamente conveniente de expre¬ 
sar las fluctuaciones de una membrana, lo de curvatura 
media k de una membrana. En los diferentes métodos 
se aprovecha el hecho de que valores diferentes de k ha¬ 
cen que la membrana sea más o menos flexible y por lo 
tanto tenga diferente susceptibilidad a las fluctuaciones 
térmicas del solvente. Así, una membrana muy flexible 
fluctúa más que una membrana rígida, lo cual ocasiona 
consecuencias medidles en varios fenómenos físicos. Por 
ejemplo, la dispersión de radiación (luz, rayos X o neu¬ 
trones) es diferente si la membrana es flexible o si es 
rígida. De igual forma, cuando una vesícula gigante es 
aspirada con una micropipeta, el cambio aparente en su 
área es diferente, dependiendo de si la membrana que 
la forma es flexible o no. Analizaremos este último ex¬ 
perimento con más detalle en la siguiente sección para 
mostrar una manera de medir k. Para ello, primera¬ 
mente derivamos algunas expresiones matemáticas para 
las fluctuaciones de una membrana. 


D. Energía de Helfrich en la parametrización de 
Monge 

Por simplicidad, vamos a analizar el caso de una mem¬ 
brana que en promedio se encuentra horizontal respecto 
a un plano de referencia pero que presenta pequeñas fluc¬ 
tuaciones térmicas en la dirección perpendicular al plano. 
La posición sobre el plano de referencia se describe con un 
vector de posición bidimensional r = (x, y). A un tiempo 
dado, la conformación de la bicapa puede caracterizarse 
mediante la altura en cada punto z = h{r) (Ver Figura[7|). 
A la descripción matemática de la membrana realizada 
así, se le llama “parametrización de Monge” y es una de 
las maneras más simples de parametrizar una superficie. 
Supondremos que la membrana tiene dimensiones L x L 
en el plano de referencia. Además, por simplicidad, con¬ 
sideremos una bicapa con dos monocapas idénticas de 



fórmula: 



Figura 7. Representación esquemática de la parametrización 
de una membrana con pequeñas fluctuaciones térmicas. 


modo que cq = 0. 

Lo primero que veremos es la forma que adopta la cur¬ 
vatura media en la representación de Monge. Comence¬ 
mos por analizar un caso más simple: una curva uni¬ 
dimensional en un espacio bidimensional. Dicha curva 
puede pensarse como la trayectoria de un móvil, R{x), 
donde el parámetro x puede pensarse momentánemente 
como un tiempo. La “velocidad” del móvil sería v = 
R'{x) y su “aceleración”, a = R"{x). El apóstrofe denota 
derivada con respecto a x. La aceleración centrípeta es 
igual al cuadrado de la velocidad por la curvatura en el 

punto en cuestión, Oc = ^ = Por tanto, para 

calcular la curvatura necesitamos calcular de manera in¬ 
dependiente la aceleración centrípeta. Esto lo haremos 
aprovechando que, a su vez, Oc es igual a la componente 
de la aceleración perpendicular a la velocidad. Esta com¬ 
ponente puede ser obtenida recordando que el área A del 
paralelogramo formado por los vectores R' y R" es igual 
al producto de su base, ||i?'(a:)|| por su altura, Oc, de 
donde Oc = Ahora recordemos que el área (ori- 

entada) del paralelogramo formado por los vectores R' {x) 


y R"{x) es igual al determinante det 


R'{x),R"{x) 


(el 


signo del determinante está relacionado a la orientación 
de un vector respecto al otro). Por lo tanto, la acel¬ 
eración centrípeta está dada por Oc = ^ 

Despejando, la curvatura en dicho punto es entonces: 
c = — = _ Esta fórmula es general para 

c ||it (íC)ll 

cualquier curva. 

Para hacer la analogía con el caso que queremos es¬ 
tudiar podemos pensar en una curva que es casi recta, 
excepto por pequeñas fluctuaciones, y que por lo tanto 
puede representarse por el análogo de la parametrización 
de Monge, R{x) = {x,h{x)). Es ahora claro que, a 
primer orden en h y sus derivadas, la curvatura es: 


„ _ det[(l,fe,,),(0,/txrr)] L 

Análogamente, en el caso tridimensional que nos in¬ 
teresa, la curvatura media de una superficie está dada, 
a segundo orden en h y sus derivadas, por la siguiente 




donde representa el operador laplaciano, el cual cor¬ 
responde a la segunda derivada del párrafo anterior. Por 
su parte, el elemento diferencial de área está dado por 


+ hl + hydxdy 


1 + i (Vhf 


dxdy. 


Sustituyendo 


estas expresiones en la ecuación (¡3]), tenemos que la en¬ 
ergía de Helfrich para membranas débilmente curvadas 
puede escribirse en forma aproximada como: 


+ aiVhf'^dxdy, 


( 6 ) 


donde estamos ignorando el término constante. Esta 
ecuación supone además que las fluctuaciones térmicas 
afectan la forma de la membrana pero no provocan cam¬ 
bios en su topología, por lo cual se ha ignorado el término 
gaussiano, que como vimos permanece constante mien¬ 
tras no cambie la topología. Esta manera de escribir 
la energía elástica de Helfrich es matemáticamente con¬ 
veniente para describir las fluctuaciones de forma en la 
membrana, lo cual ha llevado a concebir métodos exper¬ 
imentales para medir la constante elástica k. 


E. Fluctuaciones de forma de una membrana 

Para estudiar las fluctuaciones de una membrana, va¬ 
mos a considerar sus diversos modos de deformación. 
Para ello, haremos una expansión en serie de Fourier de 
la función que describe la forma de la membrana en la 
parametrización de Monge, h(r). Esta expansión queda 
de la siguiente forma: 

h{r) ='^hf^-exp{iq-r} (7) 

donde q = ^ {nx,ny) y ny son números enteros. 
Los coeficientes pueden verse como “coordenadas” en 
el espacio de Fourier y representan la contribución a las 
fluctuaciones de cada modo de deformación. 

Para utilizar esta expansión en la expresión de la en¬ 
ergía de Helfrich (ecuación [ü]), derivamos la expresión an¬ 
terior para obtener el gradiente y el laplaciano de h{f): 

{yhf = - ^ (qi • (fa) exp{z(gi -b ga) • r\, (8) 

^ qlqlh^^h¡¡^ exp {i{qi + ga) • r}. (9) 

Con estas ecuaciones, la energía libre de Helfrich, m, 
puede reescribirse como 
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- a {qi- q 2 )] exp {i{qi + ¿ 12 ) ■ r}(fr. 


( 10 ) 


Esta última expresión puede escribirse en una forma más 
simple intercambiando la integral con la sumatoria y us¬ 
ando la representación de Fourier de la delta de Kro- 
necker: 


/ exp{—iq-v}d^v = L^6^fi. ( 11 ) 

Jlxl 

¿2 

De aquí, F = — X)ij ^ 9^-9 (^9^ + ■ Por otra parte, 

el hecho de que h{f) sea real y la ecuación ([T]) nos llevan a 
que los modos complejos de Fourier son tales que h_^ = 
ht. Todo esto finalmente nos lleva a que la energía de 
Helfrich es: 


F {Kq'^ Faq^). (12) 

9 

Observemos ahora que la energía obtenida es armónica 
en estas nuevas “coordenadas” \hq\. Si ahora se aplica 
el teorema de equipartición de la energía (Ver Apéndice 
El, las amplitudes cuadráticas medias de los diferentes 
modos de oscilación en la membrana están dadas por: 




KbT 

(^(74 -p ( jq ^)' 


(13) 


Comúnmente a esta función se le da el nombre de es¬ 
pectro de fluctuaciones o factor de estructura. Vemos en 
esta ecuación que, como era de esperarse, la amplitud de 
las fluctuaciones de la membrana aumenta al aumentar la 
temperatura. Igualmente, dado que la amplitud de cada 
modo es inversamente proporcional al módulo elástico de 
curvatura media, k, vemos que las fluctuaciones son más 
grandes mientras más pequeño sea este módulo, y vice¬ 
versa. 

Por otra parte, dos modos diferentes están desacopla¬ 
dos y, por tanto (Ver Apéndice El, 


(/ig-,/ig-,) = 0, (14) 

si qi ^ —ú* 2 - Utilizando este hecho es posible ver que 
la amplitud total de la fluctuación de la membrana en 
el espacio real es simplemente la suma sobre todos los 
modos individuales, es decir, 


F) = Y.F¿>- 
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Esta relación es muy importante pues su evaluación per¬ 
mite tener una idea sobre los parámetros de los que de¬ 
pende la amplitud de las fluctuaciones de forma en la 
membrana. Para calcular la amplitud total en forma 


aproximada la suma se puede transformar en una inte¬ 
gral: 






2 

qdq. 


Cabe hacer notar que en esta ecuación, los límites de in¬ 
tegración se deben elegir de tal manera que se sumen las 
fluctuaciones que tienen sentido físico. Por ejemplo, no 
tiene sentido considerar fluctuaciones cuya longitud de 
onda sea mayor que la dimensión mayor de la membrana 
(L). De igual forma, no tiene sentido considerar fluctua¬ 
ciones más pequeñas que la escala molecular. Para esta 
última escala espacial, como distancia de corte se puede 
tomar una distancia del orden del espesor de la mem¬ 
brana, pues este tamaño está determinado por la longi¬ 
tud de las moléculas anfifílicas. Con estos razonamientos, 
con una buena aproximación se pueden proponer límites 
de integración dados por: qmax = 27r/a y qmin = 27r/L, 
donde a es del orden del grosor de la membrana y L es 
la longitud total de ésta. 

Tomando en cuenta las consideraciones expuestas, 


= 27r 



es decir, 


pQmax 

qmin 


KbT 

{nq'^ + (jq^) 


qdq, 


(F) 


KbT ql,aa, {ql^jnK- + cr) 

{q?naxl^ + ' 


(15) 


En particular, cuando qK^ ^ dmax y o" —>• 0, la am¬ 
plitud total toma la forma: 


(F) 


KbT 


(16) 


Vemos entonces que el cuadrado de la amplitud to¬ 
tal de las fluctuaciones es (en promedio) proporcional a 
la temperatura e inversamente proporcional al módulo 
elástico de curvatura media, k, lo cual corrobora la idea 
intuitiva de que mientras más flexible es la membrana, 
mayores son sus fluctuaciones. 


F. Medición experimental del módulo de curvatura 
media de una vesícula mediante manipulación con 
micropipetas 

Hemos dicho que existen diferentes métodos experi¬ 
mentales que permiten medir el módulo elástico de cur¬ 
vatura media, k. En esta sección comentaremos uno de 
ellos, basado en la técnica de aspiración de vesículas con 
micropipetas. 
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Este método aprovecha que las vesículas gigantes, 
preparadas con el método de electroformación, son fácil¬ 
mente observables con ayuda de un microscopio óptico 
[¡10. En general, las vesículas tienen una aparien¬ 
cia de esferas perfectas. Sin embargo, debido a que la 
membrana que las conforma es flexible, su superficie pre¬ 
senta una cierta corrugación cambiante con el tiempo 
debido a las fluctuaciones térmicas [1^. Esas pequeñas 
desviaciones respecto de la esfera se aprovechan para la 
medición experimental del módulo de curvatura media k 
de la membrana. 

Para ello se utiliza la técnica de micromanipulación 
o microaspiración con micropipetas. En este método, 
vesículas gigantes preparadas con el método de electro- 
formación (tamaño del orden de 20 mieras) se aspiran con 
un capilar muy delgado de vidrio (diámetro del orden de 
5 mieras). Mediante el capilar se aplica una presión de 
succión a la vesícula (Figura [S]). Como la vesícula es de¬ 
formadle, parte de ella penetra en la sección cilindrica 
de la micropipeta. El área de la vesícula cambia y su 
forma ya no es exactamente esférica, sino que se debe 
considerar el trozo de membrana dentro del capilar. Ev¬ 
identemente la distancia de penetración depende tanto 
de la presión aplicada como de la elasticidad de la mem¬ 
brana. Así, evaluando el cambio del área de la vesícula (a 
partir de fotografías), y conociendo la presión aplicada, 
puede determinarse el módulo de curvatura media de las 
membranas, k. 

Tanto el cambio de área de la membrana, como la 
tensión de la vesícula pueden escribirse en términos de 
parámetros controlables o medibles en el experimento. 

Por ejemplo, utilizando argumentos de equilibrio ter- 
modinámico, Henriksen e Ipsen mostraron que la 
tensión de la membrana, tr, está relacionada a la presión 
de succión, Ap, mediante la fórmula: 


a = 


ApR 


'P 


rCy 


(17) 


donde Rp es el radio de la pipeta y R^ el radio de la 
vesícula. 

Para evaluar el cambio del área de la vesícula al 
someterse a una presión de succión, recordemos que la 
membrana presenta fluctuaciones térmicas que la defor¬ 
man ligera pero constantemente. Una representación es¬ 
quemática de estas fluctuaciones en un momento dado 
se muestran en la Figura [Ü] En estas condiciones, el 
efecto de una presión de succión puede analizarse en dos 
regímenes cualitativamente diferentes, los cuales depen¬ 
den de la magnitud de la presión de succión y, por lo 
tanto, de la tensión a: 


1. Cuando la tensión ct a la cual está sometida la mem¬ 
brana es pequeña, el efecto de la presión aplicada 
es reducir la amplitud de las fluctuaciones térmicas, 
alisando la forma esférica de la vesícula, ver Figura 
m Entonces, al penetrar la vesícula en el capilar, 
el área de la membrana permanece constante. El 


aparente incremento de área (cilindro de longitud 
L en el capilar) viene de la supresión de las fluctua¬ 
ciones de curvatura alrededor de la forma esférica 
(Figura [5]). Evidentemente que este efecto está di¬ 
rectamente relacionado con la flexibilidad o rigidez 
de la membrana. Si la presión de succión está en 
este régimen es posible medir el módulo de cur¬ 
vatura media de la membrana, n. 

2. Cuando la tensión a aplicada a la membrana es 
grande, esto trae como efecto, además de desapare¬ 
cer las fluctuaciones de curvatura de la membrana, 
un estiramiento del área de la misma. En este caso, 
el área de la membrana se incrementa. También en 
este régimen es posible medir lo que se conoce como 
la “constante de compresibilidad de área”, Ka, que 
es el equivalente bidimensional a la compresibilidad 
volumétrica de un fluido ¡2^ . 

Aquí estamos interesados en mostrar cómo es posi¬ 
ble medir el módulo de curvatura media, k. Por esta 
razón concentraremos nuestra discusión en el régimen de 
pequeñas tensiones. Vamos a analizar con detalle el caso 
más sencillo de una membrana plana, aunque fluctuante, 
de dimensiones L x L. El resultado Anal es idéntico 
para una vesícula esférica pero es necesario mayor tra¬ 
bajo matemático para deducirlo. 

Primero consideremos que el área real de la membrana 
(incluyendo fluctuaciones térmicas) no es exactamente el 
producto de sus dimensiones, L x L. Esto es debido a 
que, por las fluctuaciones térmicas, la membrana no es 
perfectamente plana, sino que en cada instante presenta 
cierta rugosidad. El área real puede calcularse integrando 
el elemento de área que describimos en la sección anterior: 

A= í Jl + {Vh{r)f<fr 
JlxL ^ 

- Jl ^ + ^ [Vñ(f)]^^dV, 

donde, como en la sección anterior, h (r) es la altura de la 
membrana respecto al plano de referencia y su expansión 
de Fourier está dada por la ecuación 0. Obteniendo el 
gradiente encontramos que 


AA = -- 




’LxL 


91,92 


donde AA = A — Aq es la diferencia entre el área real 
de la membrana. A, y el área aparente al microscopio 
(es decir, el área de la membrana plana de referencia), 
Aq = Utilizando la ecuación (HU esta expresión se 
simplifica: 

q 


Podemos ahora tomar el promedio de AA y utilizar la 
ecuación m de la sección anterior para obtener el si- 
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Figura 8. Fotografías tomadas en un microscopio óptico de una vesícula de SOPC aspirada dentro de una micropipeta. La 
presión de aspiración es mayor en la foto de la derecha. Esto se refleja en la longitud diferente de la sección de membrana 
dentro de la pipeta (flechas oblicuas). Fotos cortesía de Gerardo Paredes, estudiante de doctorado en la Universidad de Sonora. 



Exceso de área 


Aa- el área aparente de la membrana correspondiente a 
la tensión a y definimos la siguiente cantidad: 

^Afj = Areal ~ Aa-, (19) 

Apliquemos ahora la ecuación (IT^ para dos casos dife¬ 
rentes: 


Figura 9. Representación esquemática de una membrana sin 
fluctuaciones térmicas y de una membrana con fluctuaciones. 
Al aspirar con una micropipeta la vesícula fluctuante, las fluc¬ 
tuaciones se suprimen. El exceso aparente de área (porción 
cilindrica de membrana dentro de la pipeta) es el área asoci¬ 
ada a las fluctuaciones de forma antes de la aspiración. 


guiente resultado: 


1. Si la porción de la vesícula no está sujeta a tensión 
externa (cr = 0 en la ecuación [TSl) tenemos que: 


AAq KbT Í óIq ] 

Aq Sttk 1 J 


( 20 ) 


2. Si la porción de la vesícula de interés está sujeta a 
una tensión ct tal que (27r/L)^ < ct/k < {27:¡aY, 
entonces tenemos que: 


9 


KbT 
-b cr 


_ KbT ¡ 

Acr Sttk \ a^a J 


( 21 ) 


Luego, siguiendo el mismo razonamiento para llegar a la 
ecuación m, tenemos que 

(AA) ^ KbT í ql,a^K + <j \ 

Aq Sttk \ qK^n + a} 

o equivalentemente, 

(AA) _ KbT / (27r/a)^ -b a/n \ 

Aq Sttk \ {27:/L)"^ + a /k ) 


Ahora, nótese primero que 

Areal — AAq -b Ag = AA^- -b A^, 
y, por lo tanto, 

iL = + i 

Ag Ag Ag 

Así, sustituyendo en esta última relación las ecuaciones 
(Hni) y (El]), se sigue que 


La ecuación dm fue deducida para una membrana 
plana. Sin embargo, es igualmente válida en el caso de la 
vesícula aspirada mediante micropipetas [1^ . Ahora 

bien, en dichos experimentos es posible observar como se 
modifica el área aparente de la vesícula inicialmente libre 
de tensiones cuando se le aplica una tensión pequeña. En 
todo el proceso, tanto antes como después de aplicar la 
tensión (siempre que esta última sea pequeña), el área 
real de la membrana, A^eaii se mantiene constante, como 
se explicó en la Figura |S| En lo que sigue denotamos por 


Aa 

Aq 


1 -b /31n 



^4^1n 




donde (3 = Luego, factorizando el término ^ y 

simplificando, obtenemos finalmente que 


Ag 


1 -b /31n 


1 -b ,51n 


l i 

a^cr J 


( 22 ) 
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En el experimento es posible medir directamente las 
áreas aparentes Aq y A^, y calcular la deformación 
a = = ^ — 1 que, en nuestro caso, en virtud 

de (1^ puede calcularse mediante 


1 + /31n 


a = 


í dTT^K 1 


- 1 . 


(23) 


En particular, para el caso en que /3 <C 1, condición 
que se cumple al considerar membranas rígidas como las 
formadas por fosfolípidos, la deformación puede aproxi¬ 
marse como 


a 


1 -I- /31n 



/31n 


d^a J 


- 1 , 


es decir, 


^Inj^l 

At:k ( J 


(24) 


Esta última ecuación permite calcular el cambio de 
área aparente cuando se succiona una vesícula con pre¬ 
siones bajas mediante la técnica de aspiración con mi- 
cropipetas. Como hemos dicho, en este experimento se 
inhiben mecánicamente las fluctuaciones, y el área de la 
vesícula correspondiente a la proyección dentro de la mi- 
cropipeta es igual al exceso inicial de área respecto a la 
forma esférica de la vesícula libre. Si se gráfica la tensión 
aplicada en función del cambio de área (ecuación [M]) , 
de la pendiente de la curva es posible extraer k. En la 
Figura [TÜ] se muestra una gráfica experimental (tomada 
de la referencia [l^) del logaritmo de ct vs. a, encontrán¬ 
dose que efectivamente hay un régimen lineal para bajas 
tensiones. El ajuste de la pendiente de dicha recta en los 
experimentos mencionados da una medición experimen¬ 
tal directa del módulo de curvatura media, k. El intervalo 
típi co p ara el fosfolípido SOPC es de k ~ 34 — AIKbT 

mm- 


III. FASES FORMADAS POR BICAPAS 

En la sección ITXl mencionamos que las moléculas an- 
fifílicas en solución acuosa forman diferentes estructuras 
debido a su naturaleza dual. Las estructuras básicas son 
las micelas y las membranas. Estas últimas pueden exi¬ 
stir en diferentes configuraciones, entre las cuales la más 
común es la fase lamelar. Sin embargo, las membranas 
también se agrupan en arreglos desordenados como la 
fase esponja. 


A. Transición entre la fase lamelar y la fase esponja 

Es un hecho experimental, verificado en muchos sis¬ 
temas de surfactantes, que las fases lamelar y esponja 



Deformación del área 


Figura 10. Gráfica de tensión-deformación para una vesícula 
de SOPC en solución de metanol/agua. Los puntos obtenidos 
de trazar el logaritmo natural de la tensión, a, contra la de¬ 
formación de área, a, tienen un comportamiento lineal en el 
régimen de baja tensión (0.001-0.5 mN/m). Gráfica tomada 
de [H. 



Figura 11. Diagrama de fase esquemático en términos de los 
parámetros xj4> Oue resume las características principales de 
las fases en sistemas de surfactante y agua. Gráfica tomada 

de [1^. 


aparecen en regiones vecinas del diagrama de fases. De 
hecho, hay un patrón que se repite en muchos sistemas: la 
aparición sucesiva de las fases micelar, lamelar y esponja 
al variar un parámetro fisicoquímico y, el cual puede ser 
la temperatura, la concentración de tensoactivo, la salin¬ 
idad, la concentración de un aditivo, etc. Usualmente, 
las tres fases observadas dependen de y como se muestra 
en la Figura[TT] Para valores pequeños de y se observa la 
fase micelar. Posteriormente aparecen las fases lamelar 
y esponja conforme aumenta y. 

Esta sucesión de fases intrigó durante mucho tiempo a 
la comunidad que trabaja en sistemas coloidales autoaso- 
ciativos. En particular, los investigadores se preguntaban 
cómo es posible que un cambio tan sutil en las condi¬ 
ciones fisicoquímicas (variación del parámetro y) lleve a 
una reorganización tan drástica de la estructura mem¬ 
branal. Uno de los grandes logros de la teoría de Helfrich 
consiste en que este fenómeno puede ser explicado en 
términos de las propiedades elásticas de membranas. A 
continuación veremos un modelo sencillo que explica la 
transición lamelar-esponja en estos sistemas. 
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Considerando la energía de una bicapa de moléculas 
anfifílicas es posible construir un modelo sencillo que ex¬ 
plique por qué una fase lamelar se transforma en una fase 
esponja al variar el parámetro x- Para esto se modela a 
la fase esponja como una superficie periódica mínima de 
Schwarz [2^, es decir, una superficie con curvatura me¬ 
dia H idénticamente cero. Aún cuando los experimentos 
de microscopía electrónica indican que la fase esponja es 
desordenada [1^ (a diferencia de la fase cúbica que es 
periódica), esta simplificación ayuda a visualizar dos car¬ 
acterísticas esenciales de la fase: 1) En cualquier parte 
de la estructura los dos radios principales de curvatura 
tienen similar magnitud pero signos opuestos. 2) La cur¬ 
vatura gaussiana tiene en cualquier parte una magnitud 
grande comparada con la curvatura media y tiene signo 
negativo. 

La primer característica permite despreciar a la cur¬ 
vatura media (de hecho, esta es idénticamente cero 
cuando la fase es una superficie periódica mínima). Bajo 
dichas condiciones, las contribuciones de la curvatura me¬ 
dia a la energía libre de Helfrich para la fase lamelar y 
para la esponja son básicamente iguales. Por lo tanto, la 
principal diferencia entre estas fases es debida al término 
gaussiano, como se observa en la ecuación (|4]). 

Por otra parte, es posible deducir a partir del teorema 
de Gauss-Bonnet que la diferencia en la energía de Hel¬ 
frich por unidad de área, ecuación O, entre las fases 
lamelar y esponja es del orden de Eiameiar - Eesponja = 
—ApnR, como se sigue de la ecuación ([S]). En base a esto 
podemos ver qué valores del módulo elástico gaussiano 
estabilizan la fase lamelar y cuáles la fase esponja. El 
criterio a seguir es la minimización de la energía elástica. 
Así, se espera que el valor de R que estabiliza la fase 
lamelar sea negativo, pues en ese caso, esta fase tendría 
menos energía elástica que la esponja. En cambio, para 
estabilizar la fase esponja, R deberá ser positiva. 

Para entender la transición lamelar-esponja, vamos a 
describir un modelo que predice cómo el parámetro fisi- 
coquímico X afecta el valor del módulo elástico gaussiano 
R. Para ello vamos a analizar la energía libre de Helfrich 
de una bicapa partiendo de la energía de cada una de las 
monocapas que la constituyen. 

Consideramos que las monocapas no son simétricas re¬ 
specto a la mitad de su superficie, es decir, presentan 
curvatura espontánea cq. Así, la energía de curvatura 
para cada monocapa puede escribirse como: 



Figura 12. Estructura esquemática de una bicapa. La su¬ 
perficie incompresible de cada monocapa permanece a una 
distancia e/2 de la mitad de la superficie S. 


Cuando las monocapas están unidas frente a frente for¬ 
mando la bicapa, cada una está en un estado que no 
corresponde a su curvatura de equilibrio. La presencia 
de la otra monocapa hace que necesariamente alguna de 
las dos no se pueda plegar en la dirección requerida por 
la curvatura espontánea. Para una monocapa la cur¬ 
vatura espontánea es l/i?i = cq, mientras que para la 
otra es —l/i ?2 = —cq (ver Figura [T^). Los radios de 
curvatura que tiene cada monocapa definen su energía 
elástica según la ecuación (US]). La energía libre de la 
bicapa puede ser formulada sumando las contribuciones 
de las monocapas. Así, sustituyendo la curvatura media 
en función de los radios principales de curvatura de cada 
monocapa a una distancia e/2, medida desde la mitad de 
la superficie ('Figura IT^ en la ecuación (1^ . tenemos que 




1 1 
i?l -|- ^ i?2 + ^ 



1 

2 




2 


+ Cq 


fm = [H - Cof + RmK. (25) 

Puesto que en el proceso de pasar de una bicapa plana a 
una curvada no se involucra cambio de topología, es posi¬ 
ble abandonar el término gaussiano al hacer los cálculos 
y reincorporarlo en la expresión final. Para evitar con¬ 
fusiones en este cálculo, usamos los subíndices m y b 
señalando que estamos hablando de una monocapa o una 
bicapa, respectivamente. Por ejemplo, es el módulo 
de curvatura media de una monocapa y el módulo de 
curvatura media de una bicapa. 


Haciendo un desarrollo en serie de Taylor a segundo orden 
alrededor de e/2 tenemos 


fb 




1 

lil 


+ - 


1 

R¡ 



1 



1 



Desarrollando los paréntesis y despreciando los términos 
de tercer orden y superior en los radios principales de 
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curvatura 


/b ~ 2 


2Rl 


2Rl 


R1R2 




lo que puede reescribirse como 

2 


fb = 


— —V 

Ri^ R2¡ 




coe 


1 

+ 2^" 






Considerando ahora las definiciones de las curvaturas me¬ 
dia y gaussiana en función de los radios principales de 
curvatura, tenemos que 


fb — c2^‘2KYn)H -\- 2ÍI CCQKm KcCQbvrn “t” CgAím 


Puesto que c§Krn es una constante independiente de la 
forma que adquiere la membrana y 2eH^coKm es despre¬ 
ciable en comparación con H^Km, es posible despreciar 
estos términos, con lo que se obtiene 

fb = ^{2Km)H'^ + {-ecoKm) K. (26) 


Comparando el término que involucra la curvatura gaus¬ 
siana en esta expresión con el de la ecuación ®, se puede 
concluir que 


Kb = -ecoKm- (27) 

Esta ecuación es importante pues relaciona una con¬ 
stante elástica de la membrana (Sb) con propiedades 
de las monocapas que la constituyen: su curvatura 
espontánea (cq) y su módulo elástico de curvatura me¬ 
dia (Km)- 

Esto permite entender la transición lamelar-esponja 
de la siguiente manera. Al variar las condiciones fisi¬ 
coquímicas (parámetro y) en la dirección apropiada, las 
propiedades de las monocapas cambian en consecuencia, 
por ejemplo, su curvatura espontánea. Esto ocasiona que 
el módulo de curvatura gaussiano de la membrana cambie 
(ecuación[27]), lo que, si el cambio tiene los valores apropi¬ 
ados, se refleja en un cambio de topología de membranas 
abiertas (fase lamelar) a membranas en forma de silla de 
montar a caballo (fase esponja). 

Usando la ecuación (EZH podemos ver que, en la aprox¬ 
imación en que la fase esponja es modelada como una 
superficie mínima de curvatura media nula, la energía 
libre de curvatura tiene la forma 

/ = -ecoHmK 

y por el teorema de Gauss-Bonnet, ecuación ([5]), la en¬ 
ergía sobre toda la membrana está dada por 


En esta ecuación vemos que el género de la superficie, 
p, es también un factor que influye en la energía elástica 
de la membrana. 

En esta expresión la única constante que puede cam¬ 
biar de signo es la curvatura espontánea de la monocapa. 
De donde es claro que si cq < 0 la energía mínima cor¬ 
responde a una superficie cuyo género es mayor que uno 
(pues p mayor que uno resulta en energías negativas). 
Es decir, la energía elástica favorece la fase esponja. Sin 
embargo, cuando cq > 0 la energía mínima corresponde 
a una superficie sin agujeros (p = 0) y como consecuen¬ 
cia la fase lamelar es mayormente favorecida. Vemos pues 
que este modelo sencillo permite entender, al menos cual- 
itativamene, cuándo la fase lamelar o la fase esponja son 
estables en términos de la energía elástica de las mem¬ 
branas. 


B. Interacción estérica de Helfrich 

En esta sección vamos a comentar brevemente una con¬ 
secuencia más de la elasticidad de las membranas. Como 
veremos, este efecto tiene que ver con las interacciones 
que estabilizan arreglos de membranas como la fase lame¬ 
lar o la fase esponja. 

Para ejemplificar, consideremos una fase lamelar, la 
cual presenta orden a una escala mucho mayor que el 
espesor de las membranas. La distancia media entre bi- 
capas en esta fase está determinada por las interacciones 
entre ellas. Por ejemplo, las fuerzas de van der Waals, 
que son siempre atractivas y de corto alcance, tienden a 
hacer que las membranas se adhieran entre sí y que la 
fase lamelar se desestabilice y colapse. En la práctica 
esto no sucede debido al efecto estabilizador de fuerzas 
repulsivas, como en el caso de las membranas cargadas, 
donde la interacción repulsiva es electrostática. 

Sin embargo, veremos que en sistemas sin carga tam¬ 
bién es posible estabilizar una fase lamelar con distan¬ 
cias entre membranas muy grandes, del orden de 1000 Á. 
En un principio no se entendía cual era la fuerza repul¬ 
siva que lograba equilibrar a las interacciones de van der 
Waals en este caso de ausencia de fuerzas electrostáticas. 
Ahora sabemos que esta fuerza tiene su origen en las fluc¬ 
tuaciones de forma de la membrana. 

Cualitativamente, este efecto repulsivo puede enten¬ 
derse de la siguiente manera. Cuando las membranas 
de un sistema son muy flexibles, las fluctuaciones de 
curvatura tienen una amplitud muy grande (inversa¬ 
mente proporcional al módulo de curvatura media, que es 
pequeño). Sin embargo, el confinamiento de cada mem¬ 
brana por sus vecinas en la fase lamelar reduce efectiva¬ 
mente el espacio configuracional que cada bicapa puede 
explorar con sus fluctuaciones, tendiendo a reducir la en¬ 
tropía del sistema. Es precisamente esta reducción de 
entropía, efecto no favorecido en términos de la energía 
libre, la que da lugar a la fuerza repulsiva conocida como 
interacción estérica de Helfrich. El sistema alcanza un 
equilibrio, donde la energía libre es mínima. 


F = -47recoKm(l -p). 
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Figura 13. Vista esquemática de una fase lamelar formada 
por arreglos periódicos de membranas surfactantes. 


El efecto directo del confinamiento es mejor modelado 
como una constricción a las posibles conformaciones de la 
membrana. Si la distancia promedio entre las membranas 
es d, (Figura un]), de acuerdo con la ecuación (fTBl) . el 
confinamiento será sentido por una membrana cuando su 
tamaño sea del orden de Ld, tal que 


KbT 

K 


f 2 




Podemos ahora pensar en la bicapa dividida en bloques 
de tamaño lateral Ld x Ld y suponer que cada uno de 
los bloques no es afectado por las interacciones entre las 
membranas vecinas. Así, el confinamiento no restringe 
las posibles conformaciones internas dentro de un bloque. 
Sin embargo, sí restringe el número de posibles confi¬ 
guraciones de los centros de masas de todos los bloques 
en la bicapa; esto es, restringe los posibles valores de 
Nd = {L/Ldf coordenadas configuracionales. Hagamos 
ahora la analogía con un gas ideal tridimensional de N 
partículas, cuya entropía es proporcional a 3NKb, donde 
ATb es la constante de Boltzmann. Si un número Nd de 
las coordenadas configuracionales fueran restringidas, la 
pérdida de entropía será proporcional a —NdKB- Apli¬ 
cando el mismo criterio a nuestro sistema encontramos 
que la energía libre debida a la pérdida de entropía por 
el confinamiento es 

F = -TS^KBTNd = KBT(^^^ ^^(^KbT^ . 


Por consiguiente, la energía libre de las interacciones 
estéricas por unidad de área de una membrana en la fase 
lamelar puede escribirse como 

E _ 1 (KbTV 
A K\d) 


Otras derivaciones de la fórmula obtenida pueden encon¬ 
trarse en la literatura [13, . Esta energía corresponde 

a una interacción repulsiva entre bicapas, que varía como 
el inverso del cuadrado de la distancia de separación en¬ 
tre ellas. Cuando ésta es combinada con otro tipo de 
interacción como la de van der Waals o la electrostática, 
estas interacciones determinan la periodicidad de la fase 
lamelar. Este resultado fue derivado en 1978 por Helfrich 


27¡ . mostrando que las propiedades de la fase lamelar 


son determinadas en buena medida por las propiedades 
elásticas de las membranas. 


IV. CONCLUSIONES 

En este artículo presentamos una derivación de la 
energía de Helfrich, la cual describe las propiedades 
elásticas de membranas. Según este modelo, toda mem¬ 
brana tiene dos módulos elásticos: el de curvatura media 
y el de curvatura gaussiana. El primero de ellos está 
relacionado con la rigidez de la membrana y determina 
que tanto se deforma, sin cambiar de topología, debido 
a las fluctuaciones térmicas. El módulo elástico de cur¬ 
vatura gaussiana está relacionado al tipo de topología 
que adquiere una membrana. En el texto describimos un 
método para medir el módulo de curvatura media de una 
membrana, basado en la manipulación de vesículas con 
micropipetas. Hemos visto también que la elasticidad de 
las membranas da lugar a consecuencias físicas impor¬ 
tantes, como la transición entre una fase lamelar y otra 
esponja, o la existencia de una interacción estérica en¬ 
tre membranas. Las propiedades elásticas de membranas 
juegan un papel importante en la forma que adquieren 
las membranas biológicas en la célula. Diversos grupos 
de investigación han explotado estos conceptos para pre¬ 
decir con cierto éxito la forma tanto de células como los 
glóbulos rojos [2^, o de organelos como el complejo de 
Golgi [l^. Dejamos para un trabajo posterior una de¬ 
scripción de algunos de estos resultados. Cabe decir que 
el modelo presentado es relativamente simple. De hecho, 
sólo considera membranas fluidas con pequeñas fluctua¬ 
ciones, despreciando efectos como el citoesqueleto de la 
célula, las deformaciones laterales, las proteínas trans¬ 
membranales y demás complejidades presentadas en las 
membranas biológicas. Sin embargo, en la literatura hay 
extensiones del modelo, las cuales rescatan como límite 
el modelo de Helfrich [30l - l^ . Sin embargo, pese a sus 
limitantes, el modelo de Helfrich ha sido muy exitoso y es 
una referencia fundamental en el estudio de membranas 
biológicas. 
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Apéndice A: El teorema de equipartición de la 
energía 

Consideremos un sistema compuesto por un gran 
número de partículas N (que, por simplicidad, supon- 
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dremos tienen idéntica masa m). El estado microscópico 
del sistema está caracterizado por las (!>N coordenadas de 
las posiciones (ql, q 2 , y los momentos (pi,p 2 , ■■■Pn) 

de las partículas. Supongamos que el sistema está en 
contacto con un baño térmico a temperatura T con el 
cuál intercambia energía (en forma de calor) pero no 
volumen ni partículas. Uno de los teoremas fundamen¬ 
tales de la mecánica estadística clásica establece que 
la densidad de probabilidad de que el sistema se en¬ 
cuentre en un microestado particular con energía E = 
H {(fi, ...q^\p{, ..-Pn) está dada por: 

g-/3íf ( 9 I. ■ ■-w ;pl. ■ ■-piv) 

P [qi, -qN-^Pl, -pn] = J 

(Al) 

donde ¡3 = (KbT) y Kb es \a constante de Boltzmann. 
Esta es la llamada “distribución canónica”. La energía 
cinética de la i-ésima partícula es pf/ (2m), por lo que 
es fácil ver de la ecuación (lAll) que la energía cinética 
promedio de cada partícula es: 


/^\ ^ I d^pjpí/ i^rn)]e 2m ^3 

\2m/ ^ ^ 2 ^ 

f cEpe 2m 

De aquí se sigue que la energía promedio correspondiente 
a cada coordenada de momento es [KbT) /2. Supong¬ 
amos ahora que nuestro sistema particular consiste en un 
conjunto de osciladores armónicos desacoplados, de man¬ 
era que la energía correspondiente a cada coordenada (o 
“grado de libertad”) tiene la forma: fca:f/2, con k una 
constante. Utilizando nuevamente la ecuación (ED es 
fácil ver que: 



kx1 

f dxi [kx1/2] e ^ 2 

kxf 

f dxiC _2_ 



Concluimos que cuando la energía es cuadrática en al¬ 
guna de las coordenadas (ya sea de posición o de mo¬ 
mento) obtendremos que la energía promedio correspon¬ 
diente a ese grado de libertad es igual a (KbT) /2. Este 
resultado se conoce como el Teorema de Equipartición de 
la Energía. 

Otro resultado que se sigue de la distribución canónica 
y de que las coordenadas Xi están desacopladas en la 
expresión para la energía, es que dos coordenadas difer¬ 
entes, Xi y Xj, son estadisticamente independientes. En 
efecto. 


{XiXj) 


kxf kx'j 

J dxiXiC 2 J dxjXje 2 


(A2) 


f dxiC 2 J dxjC 2 

Regresemos al caso más complejo de la membrana plana 
con pequeñas fluctuaciones térmicas de la sección llT^ De 
acuerdo a la ecuación (Un, la energía de la membrana es 
cuadrática en las coordenadas (en el espacio de Fourier) 
h^. A cada uno de los modos de Fourier (desacoplado 
a los otros modos) corresponde una energía cuadrática 
igual a -|- crg^) |/ig-p/2. En analogía con los casos 

más simples discutidos arriba, la energía promedio por 
modo de Fourier será: 


{nq^ + aq^) 


\kbT^ 


lo que nos lleva directamente a la ecuación (1131) . Del 
mismo modo, en analogía con el resultado de la ecuación 
(IA2I) , obtenemos para el caso de la membrana la ecuación 

(lil¬ 
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